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1.  VIBRATIONAL  SPECTRUM  OF  A  ONE-DIMENSIONAL  CHAIN  WITH 


RANDOMLY  DISTRIBUTED  IMPURITY  SPRINGS 

by 

Hin-Ghiu  Poon  and  Arthur  Bienenstock 
Division  of  Engineering  and  Applied  Physics 
Harvard  University,  Cambridge,  Massachusetts 

ABSTRACT 

The  frequency  spectrum  of  a  one -dimensional  lattice  containing 
randomly  distributed  impurity  springs  has  been  evahiated  to  first  order  in 
the  concentration,  a,  of  impurity  springs.  It  is  shown  that,  with  some 
mathematical  manipulation,  the  solution  can  be  placed  into  correspondence 
with  the  solution  of  Langer  for  the  analogous  problem  of  isotopic  impurities. 
The  virtual  crystal  approximation  which  yields  the  correct  elastic  constants 
requires  an  effective  spring  constant,  7,  which  is  given,  in  terms  of  normal 
spring  constant,  7,  and  the  impurity  spring  constant,  7',  by  the  relation 

i/t  =  [U-q)/?]  +  q/V  • 


I.  INTRODUCTION 

In  recent  years,  considerable  progress  has  been  made  in  the 
calculation  of  the  vibrational  properties  of  solids  with  isotopic  impurities  [l]. 
For  a  variety  of  reasons,  however,  there  has  been  less  consideration  of 
problems  in  which  the  interatomic  force  constants  are  varied. 
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Such  problems  appear,  at  first,  to  be  unrelated  to  actual  physical 
situations.  In  the  disordering  of  two  alloys  of  the  beta-brass  type,  CuZn  and 
CoFe,  however,  the  constituents  have  almost  identical  atomic  masses.  The 
primary  result  of  disordering  appears  to  be  a  change  of  force  constants.  Since 
the  change  of  the  vibrational  spectrum  wi4h  order  influences  the  equilibrium 
state  of  order,  it  will  be  necessary  to  understand  these  changes  before  a 
complete  understanding  of  the  order-disorder  process  is  possible. 

First  steps  in  this  direction  have  been  made  by  Wojtowicz  and 
Kirkwood  [2]  as  well  as  Oguchi  and  Hiroike  [3].  In  both  cases,  however,  there 
are  approximations  which  considerably  limit  the  applicability  of  the  results. 
While  it  is  highly  likely  that  any  solution  of  these  problems  will  be  approximate, 
there  is  reason  to  believe  that  more  satisfactory  approaches  may  be  found. 

For  that  reason,  we  have  turned  our  attention  to  the  formal  solutions  of  certain 
simplified  problems.  It  is  hoped  that  these  solutions  will  lead  to  some  insight 
that  will,  in  turn,  lead  to  better  approximations.  This  first  solution  of  a  one¬ 
dimensional  problem  has  been  solved  to  develop  useful  mathematical  procedure. 
Solutions  of  problems  of  higher  dimensionality  have  been  obtained.  They  will 
be  presented  when  they  are  understood. 

The  problem  which  is  treated  here  is  that  of  a  linear  chain  of  atoms 
v/ith  the  same  mass  and  nearest  neighbor  interactions.  Randomly  positioned 
springs  are  replaced  by  springs  which  have  a  different  force  constant,  and 
the  spectral  distribution  function  in  calculated  to  first  order  in  the  concentration 
of  the  impurity  springs.  The  mathematical  techniques  used  are  almost  identical 
to  those  of  Langer  [4]  who  obtained  an  exact  solution,  to  first  order  in  the 
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concentration,  for  the  spectral  density  function  of  a  one-dimensional  chain 
of  identical  a*cms  in*o  which  isotopic  impurities  are  randomly  placed.  Indeed, 
we  proceed  by  showing  that  the  problem  considered  here  can  be  placed  in  the 
same  form  as  the  problem  which  Langer  considered,  and  then  use  his  results. 

The  corresponding  order-disorder  problem  is  both  simple  and 
uninteresting.  Consider  a  chain  consisting  of  two  distinct  chemical  species, 

A  and  B,  of  the  same  mass,  which  alternate  in  position  in  the  ordered  ground 
state.  Let  the  ordering  process  and  the  vibrational  spectrum  be  controlled  by 
nearest  neighbor  interactions,  such  that  the  spring  constant;  between  A-B  pairs 
is  different  from  that  between  A -A  and  B-B  nearest  neighbor  pairs.  Let  the 
force  constants  between  A-A  and  B-B  pairs  be  identical.  Such  a  situation 
would  be  obtained  from  the  polar  model  of  Mott  [5]  or  Harrison  and  Paskin  [6] 
if  their  interactions  were  of  sufficiently  short  range  (which  ‘hey  are  not)  and 
if  the  rep^Bive  interactions  were  all  the  same.  The  first  configurational 
excited  state  of  such  a  system  is  one  in  which  there  is  an  A-A  or  B-B  pair 
somewhere  along  the  chain,  but  perfect  short  range  order  elsewhere.  That 
is,  there  is  a  simple  antiphase  boundary.  The  corresponding  vibxational 
spectrum  would  be  one  for  an  isolated  single  impurity  spring  constant  placed 
at  random  in  the  chain.  The  next  excited  state  which  maintained  the 
composition  would  consist  of  two  such  antiphase  boundaries.  One  would  be 
between  an  A-A  pair  and  another  between  a  B-B  pair,  so  that  they  would 
be  separated  from  one  another  by  more  than  a  nearest  reighbor  distance. 

(This  restriction  can  be  lifted  by  the  choice  of  suitably  chosen  grand 
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canonical  ensemble. )  It  follows  that  the  randomly  placed  impurity  spring 
constant  problem  can  be  brought  into  close  analogy  with  the  one-dimensicnal 
order-disorder  problem  where  the  concentration  of  impurity  springs  is 
analogous  to  the  decrease  of  the  short  range  order  parameter  *om  unity. 

It  should  be  noted  that  the  choice  of  randomly  positioned  impurity  springs 
makes  this  solution  applicable  to  the  order-disorder  problem,  rather  than  the 
impurity  problem.  That  is,  one  might  be  concerned  with  the  effect  of  randomly 
distributed  impurity  atoms  on  the  vibrational  spectrum.  The  major  feature  of 
these  impurities  might  very  well  be  the  changing  of  force  constants,  rather 
than  the  changing  of  the  mass.  To  treat  this  problem,  however,  one  would  be 
obliged  to  consider  randomly  distributed  pairs  of  impurity  springs  on  the  sides 
of  the  impurity  atom,  rather  than  randomly  distributed  single  springs. 

II.  FORMULATION  OF  THE  PROBLEM 

The  formulation  of  the  problem  presented  here  is  similar  to,  and  whenever 
possible  and  practical,  identical  to  that  presented  by  Langer  [4],  We  consider 
a  very  long  chain  in  which  every  atom  has  the  same  mass.  There  are  two  types 
of  force  constants  present  in  the  chain,  y  and  y '  .  The  concentration  of 
impurity  force  constants,  y '  ,  is  q  ,  where  q  «  1  The  unperturbed  chain 
has  q  =  0. 

For  a  particular  distribution  of  springs,  y  and  y'  along  the  chain,  the 
equations  of  motion  are 

mV y t,i-i  +  r»>l+1  (x/+i  -x/>  • 


a) 
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Here  x^  is  the  displacement  of  the  atom  at  the  site  I  ,  and  ^  ^  is  the 
force  constant  acting  between  atoms  at  site  f  and  1-1.  The  equations  of 
motion  can  be  rewritten  as 

1  ‘  X(>  -1'(X1+1  -Xl> 

=  (•>.,,  ,_1  -  T)  <x,.i  •  Xi>  +  !?1,  1+1  ‘I'  )(xl+l  -Xl>'  <2 

The  transformation  to  the  normal  coordinates,  ,  of  the  unperturbed 
lattice  is 


N/2 


J7T  i,  Qk  exp 
(Nm)  k-N/2 


(3) 


where  N  is  the  number  of  atoms  in  the  chain.  Substitution  of  Eq.  (3)  into 
Eq.  (2)  yields 


(“k2-“2»«k  =  -  I  V  Qk"  <4> 

k' 

where 

=  um  lsin  k/N) |  ,  wm  =  (4y/m)1/2  (5) 

and 


*k,  k' 


(-21)  sin  Qrk'/N)  V 
Rm  L  [_j 

l 


exp  [2iri  (k'-k)f/N] 


x  ^7l,f+rT)exP  (JTlk'/N)-^  |>pl 


-  y )  exp  (-ir  ik'/N)]  . 


(6) 
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Equation  (4)  is  identical  in  form  to  Langer'?  Eq.  (4).  The  entire 
difference  between  the  two  problems  is  contained  in  the  explicit  form  of  the 
f  as  expressed  in  Eq.  (6). 

For  a  particular  configure,.  :n  of  springs,  we  can  write 


[D-V)]k(k, 


(uk 


6k,k'  + 


k,  k* 


such  that  L'q.  (4)  becomes 


=  0  . 


(7) 


(8) 


k' 

Following  Langer,  the  spectral  distribution  function  for  this  configuration, 
defined  by 


g(  w) 


lim  2u 
N-*-oo  TT 


(9) 


n 

where  the  are  the  exact  eigenfrequencies  for  the  particular  configuration 
of  springs,  is  given  by 


n!”  -±-  Im  Tr  D  (u2  +  U)  .  (10) 

e—0 

Txj  obtain  the  spectral  distribution  function  for  the  random  distribution 

of  springs  g(u),  we  must  average  over  all  configurations  which  are 

consistent  with  a  fixed  value  of  q.  As  shown  by  Langer,  it  is  sufficient  to 

2  2  2 

take  the  configuration  average,  D  (u  +  ie)  ,  of  D  (u  +  i  e).  B  (<j  +  ie  ) 


is  diagonal  in  k  ,  so  that 
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g (u) -  — 


€-C  k 


(11) 


III.  CONFIGURATION  AVERAGE 

Following  Langer,  we  obtain  from  Eq,  (7)  for  a  particular  configuration, 

6 


Dk,k'(“  >  =  .  . 


5  2 


"k  -  u 

which  may  be  iterated  to  yield 


uk  “u 


"  T  ^  \  k"  Dk* ,  k' (w^ 

-  #. i  * — ' 


(12) 


Dk,k>  >  = 


1 


Uk  -y 


lr  ^lr  lr  * 

V  >i  i» 

k.  k'  2  2  '  l_j  ]  2  27  2  27 

(^kj-u  )(Uk,  -y  ) 


Irt 

) 


®k  k' 

"T + 
yk  -y 


+  . . . 


(13) 


The  average  over  configurations  is  performed  at  this  point.  Since  the  y^s 

are  properties  of  the  unperturbed  lattice,  it  i3  only  necessary  to  obtain 

configurational  averages  of  the  products  ®k  k  *k  k  . *k  k1  ' 

^  1  1  ’  2  n’ 

proceed  to  obtain  these  averages.  The  first,  <$k  k'  >  may  l56  obtained  in 
the  same  manner  as  Langer  obtained  his  Eq.  (19).  We  get 

<  ul  >  =  ('2i^(ffk^-  <  y«cp  [2»i  (k1  -k)  f/N] 


“k,k« 


L, 

l 


X  [(7|  |+J- y)  exp  (iirk'/NM^  i_1- 7)exp  (-iTTk'/N)]  >  .  (14) 
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This  expression  is  evaluated  by  keeping  I  fixed  while  summing  over  configura¬ 
tions.  We  get  zero  contribution  from  any  configuration  in  which  y ^  or 

Yj  j  A  -  ,  Since  y ^  ^  =  y 1  with  frequency  q  ,  one  obtains 

^  =  “  qlV-Y),  (15) 

and 

<®k,k«  >  =  4  q(7'  -7)^exp  [2iril(k'-k)/N] 

f 

=  qu^X.  $k  k,  ,  (16) 

where 

^  *  (V/y>  -  i  .  (17) 

Next  consider  the  second-order  term, 

<\k1®k1,kl>  =  (”4/N^m^)  sin  (irkj/N)  sin  (jrk'/N) 

x  <£  {exp27Ti  (kj-k)  l^/N  exp  2tfi  (k' ~k, )  l^/N 

ll»i2  (18) 

x  ii  +  i-7)exp(irik1/N)  .  (7i  ^  { -y)  exp  (-irikj/N)] 

x  to'#  #  +1-7)exp(7rik'/N)  -  {y.  .  ,  -7)  exp  (-irik'/N)]}  >  . 

Terms  which  ?re  first  order  in  q  appear  in  the  configuration  average 
of  the  four  terms  represented  by  the  product  of  the  square  brackets  whenever 
the  two  springs  of  one  of  the  four  terms  are  identical.  That  is,  when  l^=  , 
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=  -  1,  or  +  1 .  This  complication  distinguishes  the  problem 

considered  here  from  that  of  Langer,  in  which  only  the  terms  contribute 

to  first  order  in  q.  This  extra  complication  forces  us  to  use  a  slightly 
different  approach  in  the  evaluation  of  the  general  product  of  the  fts.  For  the 
second-order  term,  however,  it  is  easy  to  sum  the  appropriate  terms. 

For  terms  with  ,  we  obtain 

(~4/Nm2)  sinOrkj/N)  sin  (irk'/N)  6  k  k,  qiy'-y)2, 

x  [exp  JTi  (kj+k')/N  +  exp  -  TTi  (kj+k')/N’  .  (19) 

Terms  with  =  1^-1  yield 

(4/Nm2)  sin  (irkj/N)  sin  (ffk'/N)  k,  q  (y1 -7)^  exp  [iri(k1 -k')/N]  ,  (20) 

while  terms  with  =  1^  +  1  yield 

(4/Nm2:  sin (irkj/N)  sin  (iTk'/N)  «k  k,q(Y-y)2  exp  [ -iri(kj -k')/N].  (21) 

Summing  the  contributions  from  (19),  (2C),  and  (21),  we  obtain 


<  ft 


k,ki 


^  k.  >  =  KjUk'/N)qK  6k 


0(q2) 


(22) 


It  should  be  noted  that  although  the  intermediate  steps  are  different, 
the  result  is  identical  to  Langer's  Eq.  (22). 

We  now  proceed  to  evaluate  the  configuration  average  of  the  general 


product, 
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<  , 
k>kl  k1’^2 


®.  ,  > 
k  , ,  k 
n- 1’  n 


=  (-**> 
N11  m1 


M4H4) . -(V)  -C4-) 


<  I  I . Iexp2f[(k1-k)«1+(k2-k1M2+ . t^.,1  ln  ] 


I1  ll 


x  ^Tf1>l1+l*^cxp^irlkl^"^|  |i_1-T)exp(-irik1 


/N)] 


X  ^Tl  ,  I  +r?)  exp  ^ik2/N^  ■  |  _1-T)  exp  (-irik2/N)] 

2*  2  2*  2 


x  [(7f  |  -T)  exp  (?rik  /N)  -  (y  -7)  exp  (-irik  /N)]  > 

n’  n+1  n  V  n-1  n 

where  kR  =  k'  .  Equation  (23)  can  be  written  in  the  form 


^k  ,  k  > 
n-1  ’  n 


_  ( -2i)n  “  .  . 

Mn  n  JT  Sln  (  “n"  ) 

N  m  j=l 


<  I  I  -Z  {exP^[(krk)il+--+<kn-kn-l)l„!  <24> 


I,  .  .  .  i  or,  =+l  cr  =+l 
1  n  1  —  n  — 


n 

TT  K  <7 1  t  4a  -7)  exp  [v  if.k  /N)]}  > 
j=l  J  j’  j  j  J  J 


x 
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Let  us  fix  our  attention  on  a  particular  ,  For  =  +  1,  Eq.  (24)  will 
yield  terms  of  order  q  only  if,  for  all  values  of  s  =  2,  3,  .  .  .  ,  n  and 

<r  =  +  1,  or  I  =  I.  +  1  and  a  =  -1.  Similarly,  for  a,  =Hl,  terms  of  order 
q  will  be  obtained  only  if  f  -  f.  and  <r  =  -l  or  t  =  - 1  and  cr  =  +  1. 

3  1  S  Si  S 


These  restrictions  are  summarized  by  the  condition  that 
1 


's  =  'l  +  I  (<rl  -  %  > 


(25) 

4 


Thus,  we  obtain  all  the  terms  of  o  *der  q  in  Eq.  (24)  by  substituting  Eq.  (25) 

into  it,  dropping  the  summations  over  all  I.  but  l.  =  fj  ,  but  retaining  the 

summations  over  the  <r.  .  The  result  is 

J 

k,  k.  k  , ,  k 
’  1  n-l?  n 


kn 


1  20  q  (y-T)" 


m  n  n 
N  m 


n 


fk. 

J]“  sin  (-"^-) 

j=l 


exp  [  ^-(kn-k)f1  ] 


-  *1 


(26) 


"V  y  itr*r.k.  n  . 

I  ■  •  *2.  ai  exp  * — n — )JT*j  exP  i)^<ri”lrj)+<rjkj3 

<r1=+l  <r  =+l  }=Z 


n  — 


The  summation  over  in  Eq.  (26)  can  be  performed  immediately 
This  leaves  the  summations  over  the  <r  variables.  The  coefficients  of  k 

s 

in  the  exponents  which  remain  are 

{^[(^-(Oj-o-^/N}  =  (iffo^/N)  for  6  =  1  , 


{i7>  [(<r1-o-g)+o's-(o-1-o-a  +  1)]}  =  (iTT  trs+1/N)  for  s  =  2,...,n~l,  (27) 

{i7T[((r1-(rn)+an]/N}  =  (iTTo^/N)  for  s=n. 
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Therefore,  Eq.  (26)  can  be  rewritten  as 

k,  k .  k  ,,k  > 

1  1  n-1*  n 


.  ...n  r  n  7Tk.“j 

=  3Tn  5  q(y'-y)n 

N  m  M=1  I  ’  n 

X  1  (X'j)  exP^(,r2ki+'r3k2+---+‘r„kn-l+<rlkn)'  • 

tr  ir  J 


<r, .  .  .  o- 
1  n 


Since  the  <r.  are  equivalent  variables  of  summation,  Eq.  (28)  can  be  rewritten 
J 


.  > 

^  kn-l-kn 


=  Nn:f)  n  «  W-y)n  ( JT  ,in  ff)  £  Tl'j  exP  ( 
N  m  '1=1  7  t->  j=l  J  \ 

J  rr  ir  J 


/iri  a.k. 


<r, ...  t 
1  n 


In  this  form,  the  sums  can  be  performed  readily.  Using  Eq.  (5),  the  average 


becomes 


<®kk  --*k  k  >=tkk  ik”"k2"k2'’ 

’  1  n-1’  n  ’  n  1  n-1 


2  /N1"1 


Equation  (30)  is  completely  equivalent  to  the  term  which  is  linear  in 


q  in  Langer's  Eq.  (25). 
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IV.  THE  VIBRATIONAL  SPECTRUM 
Because  of  the  equivalence  presented  above,  it  is  possible  to  use 
Langer's  results  directly.  He  finds  vibrational  modes  associated  with  each 
value  of  k  ,  and  with  frequency  given  by 

(31) 


and 


Here. 


+  Ak  "  1  Fk 


u.  q«(w  2-wlr2) 

\  _  k  m 

Ak  '  T 


k  '*'*  '"m  k 

-UV  (1-K  ) 


m 


2  .  2  2.1/2 

r  “k  *  “k(“m  -“fc  > 
lk  -  -j- - « - - - - - 


u. 


m 


•u^O  -O 


where 


(32) 


(33) 


K  =  \/(l  +  M  =  1  -  (y/y",  .  (34) 

In  these  expressions,  1^  is  a  measure  of  the  width  associated  with  ach 
phonon  of  wave  vector  k  due  to  the  lattice  disorder. 


In  addition,  the  density  of  states  is  given  by 


g(u)  =  ~  ^ 


2  ,  2  2.1/2 

2  2.  um  ’“'"m  1 

"m  +  - .  :  1/2 - 


(w  2-u2)  +  iKu 

m 


-1/2 


(35) 


Since  Langer  has  discussed  his  results  extensively  and  carefully, 


we  shall  only  draw  attention  to  certain  aspects  of  the  results  which  were  not, 
apparently,  of  primary  interest  to  him,  or  are  unique  to  the  problem  considered 
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It  should  be  noted,  however,  that  it  is  1/7  and  1/7'  in  this  work 
which  correspond  to  Langer's  M  and  M',  respectively.  This  correspondence 
become1’  obvious  when  the  different  definitions  of  X.  are  examined. 

In  the  long  wavelength  limit,  and  with  7'  a  7 

U2  «  4  sin  (-^)/Kl[(l-q)-i.+  -2,  ]  .  (36) 

Thus,  in  the  elastic  limit,  an  average  spring  constant  can  be  defined  by 

±  .  ihal  +  \  .  (37) 

—  7  7'  '  ' 

7  1  1 

The  frequency  spectrum  is,  in  this  limit,  identical  to  that  which  would 
be  obtained  if  all  the  springs  were  replaced  by  the  correctly  averaged  spring. 

It  is  important  to  note,  however,  that  this  correct  averaging,  to  obtain 
a  virtual  crystal  approximation,  is  of  1/7,  rather  than  7  itself. 

While  this  result  might  be  expected,  its  application,  as  an  approximation 
to  a  two  or  three  dimensional  system,  is  highly  suspect.  The  obvious  reason 
for  this  is  as  follows.  When  the  linear  chain  is  placed  in  tension,  the  atoms 
adjacent  to  the  impurity  springs  are  constrained  to  move  along  the  line  of 
tension.  The  analogous  situation  in  three  dimensions  is  one  in  which  tension 
is  applied  in  a  manner  which  corresponds  to  the  limit  of  a  pure  longitudinal 
mode  of  the  unperturbed  lattice.  When  the  impurity  springs  are  introduced, 
however,  equilibrium  may  be  reached  by  displacements  which  have  components 
normal  to  the  direction  of  the  tension.  In  this  case  2 q.  (37)  is  not  valid. 


*1 

.1 

rS 
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It  is  also  of  some  interest  to  examine  the  form  of  ,  as  given  in 
Eq.  (33).  It  is  easily  seen  that  =  0  for  =  0  or  .  Since 

is  positive,  it  must  have  a  maximum  somewhere  in  between  these  two  points. 
That  maximum  occurs  for 

i  y/2l  (38) 

2  (1-0 

for  all  possible  values  of  K2.  The  functional  form  of  is  shown  in  Fig,  1. 

The  maximum  is  quite  evident.  In  the  paragraphs  which  follow,  an  attempt 
is  made  to  explain  this  behavior. 

The  mathematical  origin  of  these  widths  is  discussed  by  Langer.  Each 
member  of  the  ensemble  is  not  periodic,  and  k  does  not  label  an  eigenvector. 
Thus,  a  measurement  which  fixes  k  will  yield  a  spread  in  frequencies. 

Although  k  becomes,  in  a  certain  sense,  a  good  label  for  the  ensemble  average, 
the  spread  in  frequencies  must  remain.  Mathematically,  we  see  the  spread 
in  the  form  of  the  average  spectral  distribution  function, 


g  (w)  = 


Zu 

7T 


lim 

N“*oo 

e~*0 


N 


Im 


Dk<“ 


+  ie) 


For  the  unperturbed  chain, 


(39) 


Dk  (u2  +  ie)  =  l/ffc^2  -u2  -  ie)  . 


(40) 


te  poles  of  Dk(u^)  occur  at  the  eigenvalues  associated  with  k  .  Langer's 
results  may  be  paraphrased  as  follows.  For  the  perturbed  problem,  there  are 
a  number  of  frequencies,  denoted  ^  ,  which  are  positions  of  the  poles  of 
each  D  (u2+  ic)  .  These  poles  lie  at  the  solutions  of 

K 
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“k2  -  + ci- <  i  >*° 

where,  from  Langer’s  work, 

*Kuk 


(41) 


Gk<“  >  -  --T^r 

Ku  \ 


(42) 


1  + 


"IT  /j  1  2 

k*  «k' 


to  first  order  in  q 


The  contribution  to  g(w)  from  any  one  of  the  frequencies  associated 
with  k  is  proportional  to  the  residue  at  the  corresponding  pole  of  (u^), 
which  is 


Res  {Dk(u2)}  =  [-2u  .+(»/»o)  G  <u2)|  ] 

u=uk,  i  ’  1 


-1 


(43) 


As  indicated  in  Fig.  2,  this  residue  is  largest  for  i  in  the  vicinity  of 

and  diminishes  rapidly  on  either  side  of  this  point  as  OG/8u)|y_  becomes 

^k,  i 

large. 

With  these  results  of  Langer,  the  mathematical  origin  of  the  form  of 
Fig.  1  is  readily  visualize'’  For  each  k  ,  ^  is  proportional  to  the  range 

in  u  over  which  the  ^  have  appreciable  residues.  Let  n^  be  the  number 
of  modes  for  which  the  residues  associated  with  the  .  are  appreciable, 
and  be  the  density  of  modes  of  the  unperturbed  system.  Then 

Fk  *  VK  “  "k  (um2-  “k2)‘/2  ■  <44> 
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Thus,  it  is  the  singularity  in  the  density  of  states  which  causes  the  vanishing 
of  rk  for  y^  =  wm  .  That  is,  although  n^  may  be  large  for  k  =  ^max, 
the  infinite  densi.y  of  states  causes  the  total  energy  spread  to  be  very  small. 


It  is  evident  that  F^  =  0  for  k  =  0,  The  form  of  the  increase  of 

with  increasing  k  ,  up  to  the  maximum  is  shown  more  explicitly,  however, 

2  2  2 

by  an  examination  of  Figs.  2  and  3.  In  these  figures,  Gk(y  )  and  (u  -y^  ) 

are  plotted  versus  y^  for  the  cases  to,  ^  «  y  ^  and  y,  ^  zz  y  ^  .  The 

Km  Km 

2  2  2  ° 
intersections  of  G^(y  )  with  (y  -y^  )  occur  at  the  values  of  y  equal  to 

2  2 
the  ^  .  The  residue  associated  with  tnese  y^  ^  is  appreciable,  however, 

only  at  those  values  of  y^  ^  ,  such  that  BG/Bu  is  small.  In  Fig.  2,  only 

the  central  y^  ^  satisfies  this  condition,  while  in  Fig.  3,  there  are  a  number 

of  terms  for  which  the  residue  is  appreciable.  Thus,  in  the  form  of  Eq.  (44), 

the  maximum  in  may  be  pictured  as  a  competition  between  the  slowly 

increasing  n^  and  which  increases  more  slowly  for  small  k  ,  but 


increases  rapidly  near  k 


max 
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FIGURE  2 


2.  THE  ELASTIC  SPECTRA  CF 


TWO  DIMENSIONAL  DISORDERED  LATTICES 

by 

Hin-Chiu  Poon  and  Ai'hur  Bienenstock 

ABSTRACT 

The  elastic  vibrational  spectra  of  perfurbed  square  lattice  systems  with 
nearest-neighbor  central  and  noncentral  interactions  have  been  derived.  The 
unperturbed  system  consists  of  masses  m ,  on  the  lattice  points  and  interacting 
with  force  constants,  a  ,  which  determines  the  resistance  to  compression,  and 
(3  ,  which  determines  the  resistance  to  shear  along  the  direction  (10).  In  one 
case,  the  perturbations  are  Nq  randomly  positioned  isotopic  impurities  of 
mass  m',  where  N  is  the  number  of  lattice  sites.  It  is  shown  that  the  elastic 
spectrum  for  this,  and  all  other  isotopic  impurity  systems,  is  completely 
determined  by  the  average  mass,  m  =  (1  -  q*  m  +  qm'.  In  the  other  case, 
corresponding  to  certain  order- disorder  situations,  the  constants  describing 
the  interactions  between  2Nq  randomly  positioned  pairs  of  nearest  neighbors 
are  replaced  by  a'  and  f3'.  To  first  order  in  q  ,  the  resulting  elastic  modes 
are  then  completely  determined  by  ar  average  a  ,  a  ,  where 

-  n  i  q  {a'-a  ) 

«  =  «  L 1  + - - — a — ■ — : - “T - TT7T  ■» 

a  +  U*' -  <*  )(2/ir  )tan  ' 

and  which  is  obtained  from  the  previous  expression  via  the  interchange  of 
a  and  .  The  appearance  of  (3  in  the  expression  for  a  implies  that  the 
virtual  crystal  approximation  fails.  It.  is  shown,  however,  that  two  different 
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forms  of  the  virtual  crystal  approximation  place  vounds  on  the  a  and  (3  ,  in 
accoraance  with  a  general  theory  of  Paul.  A  physical  interpretation  of  the 
results  is  also  presented. 

I.  INTRODUCTION 

In  a  recent  paper  [1]  (denoted,  hereafter,  as  P-B)  the  authors  obtained 
an  expression  for  the  vibrational  frequency  spectrum  of  a  one -dimensional 
chain  of  atoms  with  the  same  mass  and  nearest -neighbor  interactions.  Randomly 
positioned  springs  were  replaced  by  springs  which  had  a  different  force  constant, 
and  the  s  ^ctral  distribution  function  was  calculated  to  first  order  in  the 
concentration  of  the  impurity  springs.  The  mathematical  techniques  used  were 
a  modification  of  those  used  by  Lange r  [2]  to  solve  the  identical  problem  for 
the  same  unperturbed  system  in  v/hich  the  impurities  were  atoms  with  a  different 
mass.  While  the  motivation  for  the  previous  paper  was  an  attempt  to  obtain 
more  information  about  the  effect  of  disordering  on  the  vibrational  spectrum  of 
binary  alloys  such  as  CuZn  and  CoFe,  ir.  which  the  masses  of  the  two  constituents 
are  almost  identical,  the  one  “dimensional  nature  of  the  problem  limited  the 
value  of  the  solution.  Nevertheless,  the  mathematical  techniques  developed 
are  quite  helpful,  and  will  be  used  aere. 

In  this  paper,  two  related  problems  are  treated.  In  both  cases,  the 
unperturbed  lattice  is  a  square  lattice  array  of  atoms  with  same  mass,  and 
with  nearest-neighbor  central  and  noncentral  forces.  In  one  case,  the  two 
force  constants  describing  the  interaction  between  randomly  positioned  nearest- 
neighbor  pairs  of  atoms  are  replaced  by  a  different  pair  of  force  constants. 

In  the  other,  randomly  positioned  atoms  are  replaced  by  atoms  of  a  different  mass. 
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For  cc  nvenience,  we  call  the  problems  the  impurity  spring  and  impurity 
mass  problems.  In  both  cases,  the  elastic  constants  are  obtained  as  a  function 
of  concentration,  q,  to  first  order  in  q. 

The  most  striking  feature  of  the  results  are  that  while  the  virtual  crystal 
approximation  is  found  to  be  valid  fcr  the  impurity  mass  problem,  it  is  valid  only 
for  special  cases  in  the  impurity  spring  problem.  The  explanation  and  implications 
of  this  result  are  discussed. 


II.  THE  UNPERTURBED  CASE 

In  the  unperturbed  ■~ase,  we  are  considering  the  vibrational  spectrum  of 
a  square  lattice  a^ray  of  atoms  of  mass  m.  In  order  to  maintain  both  stability 
and  only  nearest-neighbor  interactions,  these  infractions  must  be  ncncentral. 
We  label  the  atom  at  the  origin  by  (0,0)  and  its  four  nearest  neighbors  by  (1,0), 
(0,  1),  (1,  0),  and  (0,  1).  The  equations  of  motion  of  the  i-th  atom  are  then 


m  x 


l,  a 


act 


,  it , 


x  ,  ,  =  0 

a 


(1) 


t\ax 

Here,  the  V  summation  runs  over  the  atom  £  and  its  four  nearest  neighbors 
and  the  a '  summation  runs  over  the  two  orthogonal  coordinates. 

The  force  constant  matrices  have  ihe  forms 


1 


C  (1,0) 

c  (0,  1) 


C<T.°>  *  -  (o*  |j) 
C<°-T>  =  ■  (o  a) 
2  '■*  '  P  »  (o  l) 


C  (0, 0, 


(2) 
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For  central  force  interactions,  P  =  0  and  the  lattice  shows  zero 
frequency  shear  modes  for  propagation  along  certain  principal  directions. 
Expressing  the  displacements  in  terms  of  normal  modes, 


[<)T  "  I 

k>  j 


(3) 


where  j  (=  1  or  2)  labels  the  two  orthogonally  polarized  modes,  and  7  labels 
a  displacement  component,  we  obtain  the  normal  mode  solutions, 


2 

wl,k 

_  _4_ 
m 

[a  sin^ 

k  a 

-Z--H3 

0  k  a 

Sin2-f-  ], 

2 

W2,k 

_  _4_ 
m 

[(3  sin^ 

k  a 

x  1 

—  +  « 

,  k  a 

sin  ]  , 

(4) 

with  eigenvectors 

Qj(k)  =  (  0  )  ,  Q2(k)  =  (®  )  .  (5) 

Equation  (5)  implies  that  for  any  propagation  vector  k  ,  the  two  normal 
modes,  1  and  2,  have  polarization  vectors  along  x  and  y  ,  respe  .tively. 
Thus,  one  of  the  two  labels,  j  and  7  ,  of  Qj ,  7  in  Eq.  (3)  is  redundant,  and 
will  not  be  used  hereafter. 
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T-* 

4>  -  , 


THE  IMPURITY  SPRING  PROBLEM 


A.  Formulation 

In  this  problem,  ’he  force  constants  ,  a  and  p  ,  associated  with  2N 
randomly  positioner  pairs  of  atoms,  are  replaced  by  force  constants  or'  and 
P'  .  To  simplify  the  equations  of  motion,  the  four  nearest  neighbors  of  an  atom 
which  are  separated  from  it  by  vectors  (1,0),  (1,0),  (0,  1),  and  (0,”I)  are 
denoted  by  '14  1,  7+2',,  (ir  3),  and  (14  4),  respectively.  The  equations  of 
motion  can  then  be  wri’ten  in  the  general  form 


!«  °) 

^2xf,  l  "  xi+i.  i  “  xi+2  l 

tel 

(0  Pj 

^2  x|,  2  ‘  X>+1,  2  “  Xi+2,  2/ 

/p  0  /  Z*tf  \  “  x<f3.  1  ~  x|+  4.  1 

\0  a)  \z^t 


,  2  ~  xf4  3,2"  X!+  4,  2) 


or  -  or 


a  a 


i ,  2/ 


Xf,  1  Xf.+ 1 ,  1 


Xi,  2  "  X*+l, 


x .  ,  -  x . .  _  .  1 
I,  i  1+2,  1 


xf,  2  *  X!+2,  2 / 


/P 

^ir  3 

0 

x<,  i 

‘  Xi+  3,  1  \ 

l 

0 

*  '®1,  £4-3  / 

X|,  2 

‘  Xi+  3,  2  / 

/p 

13 i,  /*  4 

0  \ 

S,1 

~  X|+4,  1  ^ 

\ 

0 

a-  O' ,  * 

/,  f-4 

X. 

i  f,  2 

Xir4,  2  , 

(6) 
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Equation  (6)  can  be  rewritten  in  terms  of  the  normal  coordinfxtes  of 
the  transformation 


xVy  =  L  V-]  expi-‘  - 

k 


(7) 


Since  all  the  force  constant  matrices  are  diagonal,  the  x  and  y  components  of 
the  displacement  are  independent.  Thus,  equations  of  motion  can  be  written 
simply  for  each  direction,  using  Eqs.  (6)  and  (7).  In  addition,  the  two  directions 
ai  indistinguishable.  Thus,  it  is  sufficient  to  study  one  direction  of  polarization 
over  a  complete  quadrant  of  the  propagation  plane.  For  the  x-components  of  the 
displacements,  we  have 

-U2m  x,  1+  a  (2 x^  j  -  xm ^  r  xi+2>  j )  +  p  (2 j  -  x|+3>  j  -  x|+4>  j ) 

=  rxi+i,  i)  +  {a  "ai,f+2)(xf,  rxi+2,  i* 

+  ‘^f,f+3^xf,  rx!+3,  l)  +  ^“^f,f+4)'xf,  l“xf+4,  1*  •  (8) 

By  substituting  Eq.  (7)  into  Eq.  (8),  and  using  Eq.  (4)  to  simplify  the 
resulting  expression,  one  obtains 


-m  u  Qi^>  exp  1  k*  +  m  T 

kQl^-^  exp  i-’  -f 

k  k 

=  (a  -ai  |+1)  ^  Qj(k)  (l-expikxa)  exp  ik-  Rjt 
k 

+  (cr -a^  j+2)  ^  Qj  { k)(l-exp-ikxa;exp  ik  Rj 

+  <e>-fwl  Qi'y(1  -expik^a)  exp  i  k  R< 
k 

+  {p-Pf  |+4 (1-exp-ik  a)  exp  ik-  Rf  .  (9) 

k 
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Ai*--r  loth  sides  of  Eq.  (9)  by  exp-ijk*-  R|  and  summing  over 

/  ,  one  oot.-iin? 

m  N  'w/  y,  •  w“',  Qj  k1/ 

A 

-  N  Q .  k  -yp  j  <k-ks;  •  R£ 

j,  k 


Aft  •  r 
manipula  tier. 


v  [q>  .  f  ,  ?  i  -expikxa)  +  {■ a-a ,  2){l-exp-i kxa) 

P  l3;  l  *xPikya>  f  ((3  ~$t  ,+4)(1-exp-ikya)]  .  (iO) 

r  --  h  nging  k  and  k!  and  performing  some  trivial  mathematical 
s ,  Eq  (10  ::an  be  written  in  the  form 


:v'“)  U  ,  k  ) 


1  ®k!k'  Ql¥' 
k 


(ID 


where 


%  V 


<m  / 


{exp-Hk-k')  • 


[  .  ‘j  •  sir.  rk:  a  /2)  exp  ii  k‘  a/2) 

i\  i1  >  x  r  ’  x 


, 


a  '  s.:r»  \k!a -2)  exp  (-ik*  a/2) 


The  sup 


■  "f  ,•  P  s:n  (k^a/2)  exp  (i  k^a/2) 

!P,  .-p  e:r.  ,kvt/2j  exp  {-i  kU/2)J  }  .  (12) 

'  J  '  A 

pi  i  ■>.  ^  is  used  to  indicate  that  polarizations  along  x 


only  a r  ■ 
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Equation  (11)  is  identical  in  form  to  Langer’s  [2]  Eq.  (4)  and  Eq.  (4)  of 
(P-B).  The  three  problems  differ  only  in  the  form  of  the  We  shall 

follow  the  approach  described  in  (P-B)  because  it  is  closer  to  the  methods  used 
here,  although  the  basic  technique  was  originally  presented,  and  is  more 
completely  described,  in  Langer's  paper. 

For  a  particular  configuration  of  springs,  we  can  write 

[DU)(u2)]^\,  =  -  u2>  k,+  *k  k,  (13) 


such  that  Eq.  (11)  becomes 

^  [D(1>(U2)]'‘k,  Qj(k')  =  0  .  (14). 

k' 

The  spectral  distribution  function  for  his  configuration,  defined  by 


-  N-*jo 


Cj 

n 


(15) 


where  the  are  the  exact  eigenf requencies  for  the  particular  configuration 
of  springs,  is  given  by  [2] 


g  M  =  ^  TT  ^  Tr  D<1>  <“2  +  1  e) 
€—0 


(16) 


In  this  expression,  use  of  symmetry  has  made  it  unnecessary  to  include 
D(2)(:j2  +  ie)  . 

To  obtain  the  spectral  distribution  function  for  the  random  distribution 
of  springs,  g  (u)  ,  we  must  ?.verage  the  distribution  function  over  all 
configurations  which  are  consistent  with  a  fixed  value  of  q  .  As  shewn  by 
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i  j  »  p 

Langer,  it  :s  sufficient  tc  take  the  configuration  average >  T)  (u)  +  ie),  of 

ju  2  _ (l )  2 

D'  '(y  -t  i  e'*.  As  n  the  one -dimensional  problem,  D'  '(w  +ie)  is  diagonal 

in  k  ,  sc  that 


giw  = 


y  D,i<v 


4w  Jim  1  T 

■■  "  y-y  “xT  Im  /  -k-j  ■. 

7T  N-  CO  N  /_J  k  ' 

c  “*■  0  k 


i  • 


(17) 


B.  Configuration  Averages 
As  in  T ,  lu i")  may  be  iterated  to  yield 


1 


nil>  2,  _ 

Dk.  k'  ^  )  '  2  T 

-  -  i*)k  •=  w 


k,_k| 


k,k' 


V 


1 

V  "k.  k' 
> 


I 


<£. 


?:d 


k,  k 


v  * 


/  . 


-l 


w2k«>; 


k1 


2, 

U)  ) 


J 


(18) 


The  average  over  configurations  is  performed  at  this  point.  Since 

the  are  prop-mnee  cf  the  unperturbed  lattice,  it  is  only  necessary  to 

ob+ain  configura'ior.al  averages  of  the  products  I 

-1  ill »  — 2  -2>-3  -n’- 

We  proceed  to  obtain  these  a\  mages.  The  first,  < 5>k  ^,>  is  easily  obtained. 


We  need 
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*<U  > 

®k,  k1 


_  (-20 
"  TsfrrT 


(k-k1) 


x  (sin(kUt/2)[(dj  i+1  -a)exp  (ik^a/2)- (a^  |+2-a)exp(-ik^a/2)] 

+  sin  (k^a/2)[01  |+3-(3)exp-ik^a/2)-(pi  |+4-£)exp(-ik^a/2)] }  >  .  (19/ 


This  expression  is  evaluated  by  keeping  I  fixed  while  summing  over 
configurations.  No  contribution  to  the  terms  linear  in  the  a’s  is  obtained 
from  any  configuration  in  which  a ^  o  %  a ^  is  equal  to  a  .  The 

analogous  situation  holds  for  the  terms  linear  in  (3.  Denoting 
Ao  =  a*  -  a 


and 

AP  =  (3'-3,  (20) 

the  configurational  average  is 


<  ®k‘k' >  =  I  exp 1  t'-y  h 

l 

x  q  [2iAasin2(k' a/2)  +  2i  A(3sin2(k'  a/2)]  .  (21) 

x  y 

Performance  of  the  summation  over  l  yields 

<  >  =  (4q/m)  [Aas.ln2(k^a/2)  +  A0  sin2(k^a/Z)]  6k  fcI  (22) 


Next,  consider  the  second-order  term 


ARPA-18 


2-11 


<*!li  > 

k,  k,  k-.k 

—  —  i  — 


.  y  «pi!k1-k>-£l«pi(k|._kI)-al 

nr  t  „  *  1  *2 


*1’  *2 


x  <  {sin  \klx*  /2'if'ar^  #  f  j  *a;exp  (iklxa/2)-(«i  ^ +2  •0’)exp{ -ik^a/2)] 

+  sin(klya/2)[;p  ii+3"P>exP  (i  kiya/2MP|  i^+4-p)exp(-iklya/2)]  } 


x{sinik' a/2)  [[a,  ,  ,  ~af)exp  (ik *  a./2)~[at  .  ,  7 -ar)exp(-ik' a/2)l 

x  '  L  r  x  l^,  l^+Z  x 

•  sin  k'a/2}[{0  .  t , ~0)exp(ik'  a/2)-{0#  #  -0)exp(-ik'  a/2)]}>  .  (23) 

>•  J  y  /2>  x2t"’  y 


As  before,  only  terms  in  which  the  a.  .  and  0..  .  differ  from  a 
’  h)  hi 

and  0  ,  respectively,  contribute  to  the  right-hand  side  of  Eq.  (23).  In 

addition,  a?  in  sT-B),  only  the  product  terms  in  which  the  a.  and 

ll’Jl 

0.  .  represent  the  same  springs  as  the  a.  .  and  0.  .  ,  respectively, 

ll»Jl  J2’J2  l2’J2 

contribute  rc  first  order  in  q.  It  should  be  noted  that  all  the  a-terms 
correspond  to  springs  which  are  parallel  to  the  direction  (1),  while  all 
0  terms  correspond  to  springs  which  are  perpendicular  to  the  direction  (1). 
Thus,  all  prcduc*  terms  which  involve  both  a-'s  and  0's  must  vanish,  to 
firsf  order  in  q  ,  for  a  product  of  “S's  of  any  order.  As  a  result  of  these 
considerations,  we  may  write 
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®kk  ^k^k'  >  =  ~T~T  7  exp[i(k  -k)-R  Jexp[i(k<-k  )*R  ] 

-1  -1’-  N  m  1  *2 

£1>  *2 

X  ^  <  i1-H'Cr)exp  ^iklxa/2} 


-  {<*|  |  +2-a)  exp  (-iklxa/2)l  sin  (klxa/2) 

x  [(a|  f  +1-or)exp(  ik^  a/2) 

-  |  +2-a)  exp  (-ik^a/2)]  sin  (k^  a/2)  > 
2^  2 


+  <(P|  t  +3-P)  exp  (ikj  a/2) 


T’  1 


"  ^|  |  +  4'^^  CXP  ^"iklya/2)i  Sin  (klya/2' 


x  t  +3_,3)  exp  <ik’  a/2) 

* 2 9  ^2  y 

"  tP#  /  +4~P)  exP  ("ik’a/2)l  siri  (k  a/2)  >  ) 

f2»  2  y  y 


(24) 


The  important  feature  of  Eq.  (24)  which  makes  this  problem  soluble 
is  that  the  two-dimensional  average  of  Eq.  (23)  has  been  split  into  two  one¬ 
dimensional  averages  which  are  mathematically  identical  to  those  which 
have  been  evaluated  in  (P-B).  This  splitting  occurs,  for  terms  which  are 
first  order  in  q  ,  for  all  order  products  of  the  $'s.  As  a  result,  using 
Eq.  (30)  of  (P-B),  we  obtain  the  general  expression 
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C,  Calculation  of  the  Elastic  Spectrum 

The  form  of  Eq.  (25)  indicates,  as  expected,  that  is  diagonal 

in  k.  As  a  result,  we  may  follow  Langer  and  write  the  diagonal  elements 

of  H(1)  as  D where 
k 


D1[1){w2)  - 


2  2 ,  77177  2. 

Up.  -  u  +  Gk  (u  ) 


(26) 


The  configuration  average  of  Eq.  (18)  may  then  be  rewritten,  using 
Eqs.  (25.  and  >26)  as 


r  - 1  ,  „ .  n 


n=  v  ■;  i)  *4) 

L  M  n~7  n 

^  N  rn 
n 


ikx*)fx  8in2(klya/2) 

L,  2  2 

wl,k  ‘  W 


*  n  .  2 
sin 


n-  1 


hzl' 

l(v 

sin2  (klya/2)  ^ 

n-  1 

1 2 1 

2  2 

wl,k  i 

(27) 
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Since  there  are  no  terms  in  which  both  A  or  and  A  (3  appear,  Eq  (27) 
is  easily  transformed  into  the  expression 


(j)  [4q(Aa)/rn]  sin2  (k^a/2) 


,  ,  4  (A a) 

i  +  - Tf - 

mN 


Aft)  Y 

N  L 


sin  (k,  a/2) 
lx 


U1  k  ’  “ 
1  > 


[4q(A£)/m]  sin  (k  a/2) 


k  yl  k  '  W 
_1  >  -1 


The  entire  problem  is  then  reduced  to  an  evaluation  of  the 


summations 


■4rS 


sin2(kxa/2) 


k  "l.k 


j  ■  i  7. 


sin2  (k  a/2) 


k  wl,k 


Because  of  the  symmetry  associated  with  this  lattice, 


I  [a  ,  p )  =  J  (p,  a)  .  ( 

Eq.  (31)  implies  that  one  need  evaluate  Eq,  (29)  only.  Transforming  the 
right-hand  side  of  Eq,  (29)  into  an  integral,  we  obtain 


IT  /aTT/a 


T  _  lim  a 
~c— 0  47T 


sin  (k  a/2)  dk  dk 

X  XV 


(4a/m)sin  (k^/2)  +  (4p/m)  sin  (kya/2)-u>2-ie 
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In  Appendix  I,  it  is  shown  that  the  value  of  this  integral  is  given  by 

I  =  vm/2iro)tan  ^(o?/P)^2  +  (im^u^/32  a$)($/a)^  ^  (33) 

under  the  conditions 
2 

mu  /4 a  «  1 

and 

mu2'/4p  «  1  .  (34) 

These  conditions  limit  our  discussion  to  the  elastic  portion  of  the 
spectrum.  It  would  be  reasonably  simple,  but  tedious,  to  obtain  higher 
order  terms  in  the  u2  power  series  expansion  of  I.  Since  interesting 
physical  results  are  obtained  from  the  lowest  order  terms,  we  shall  not 
continue  the  procedure,  but  shall  continue  with  the  determination  of  the 
elastic  portion  of  the  spectrum. 

Substitution  of  Eq.  (33)  into  Eq.  (31)  yields 

J  =  {m/27T(3}tan  ^(p/a)1^2  +  (im2u2/32op)(a/P) 1  ^2  .  (35) 


Substitution  of  Eqs .  \  3 1 )  and  (33)  into  Eq.  (28)  yields 

(4qAar/m)  sin2  (k^a/2) 

—  1  +  l2Aa/va)±an  *  (a/p)* '  c+  (imu2Aa/8«P)(p/a)^ 


(4qAP/m)  sin4"  (k^a/2) 

1  t  (2AP/7rP)tan" 1  (p/o) 1  (imu2AP/8aP)(<r/p) 1  ^2 


f 
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As  in  one-dimension  [1, 2]  we  may  associate  the  positions  of  the  poles 

( i )  2 

of  the  function  (u  )  with  the  modes  of  vibration  of  the  lattice  with  the 

randomly  distributed  impurity  springs.  The  frequencies,  for  each  k  }  and 
with  polarization  vector  parallel  to  x  ,  are  given  by  the  real  part  of  the  roots 
of  the  equation 


2  2  .  ^(1).  2,  n 

u>k  -  w  +  '  (u>  )  =  0  . 


(3?) 


(1)  2 

A  complete  solution  of  G,  (w  )  would  lead  to  a  number  of  roots,  for  fixed 
k  ,  of  Eq.  (37).  The  existence  of  these  roots  is  a  mirror  of  the  fact  that 
k>  for  any  one  member  of  the  ensemble,  is  not  a  good  label.  By  restricting 
our  calculation  of  G^(u^)  to  small  ,  and  keeping  terms  of  lowest  order 
only  in  the  resulting  small  parameters,  we  have  suppressed  the  multiplicity 
of  these  roots.  The  resulting  expression  for  the  elastic  mode  frequencies  is 
given  by 


4 a 
m 


m  L 

The  solution  of 


q  (A«  )  _ 

+  (2/jr)(Aa)tan~1  {a/p)1/Z 

q<*P>  _ 

P+(2/7T)(A(3)tan  _1((3/a)1/2 


■1 

1 


sin2(kxa/2) 


sin2(kya/2)  . 


Eq.  (37)  also  contains  an  imaginary  part. 


(3b) 


Writing 


the  solution  as 
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the  imaginary  part 


quk 


I  (Acr)2  sin2(k  a/2) 


4(«p) 


[  1+  (2Aa/7TQ!)tan  "  *  (cr/P)* //^]2 
(Ap)2  sin2  (k  a/2) 


P[1  +(2ApAP)tanJl  (p/a) 


77 r? 

J 


(40) 


corresponds  [  1  ]  to  the  reciprocal  of  the  lifetime  of  the  mode.  It  should  be 

noted  that  F.  ,  is  of  one  order  higher  in  the  small  parameters  of  Eq.  (34) 

1  >  ± 

than  is  A^  ^  .  Thus,  although  the  low  frequency  modes  are  shifted  in  frequency, 
they  are  quite  well  defined. 

Finally,  the  density  of  states,  g(u),  as  indicated  in  Eq.  (17),  and 
calculated  in  Appendix  II,  is  gi  .  sn  by  the  expression 

tom 


g(oo) 


(41) 


Here,  a  and  P  are  defined  by  r*qs.  (57)  and  (58).  These  results  will  be 
discussed  in  Section  V,  after  the  analogous  expressions  have  been  obtained 
for  the  isotopic  impurity  problem. 


IV.  THE  ISOTOPIC  IMPURITY  PROBLEM 
A.  Formulation 

In  this  problem,  the  unperturbed  lattice  is  again  the  square  lattice 
discussed  in  Section  II.  In  this  case,  the  Nq  atoms  at  randomly  positioned 
sites  are  replaced  by  isotopes  of  mass  m*.  The  mass  of  the  isotope  at 
site  i  is  denoted  by  m^.  The  equation  of  motion  of  the  atom  at  this  site 
is  given  by 
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-u 


+.l  m 

I 

r1, 2  i 

m  0  j3 
\  1 

l 

xi,  1"  xf+ 1 ,  1 

li,  2"  Xf+1,  2 


\  +  J_ 

( a  o\ 

/ 

I  m 

0  P 

/ 

^  i 

xi,  r  xi+2,  i 

xi,2'  Xf+  2,  2 


j_  /  P  o 

+  m  |  C  a 


x|+3,  l' 

+  — 

Ip  0] 

X<+3,2j 

m 

0  a 

l  i 

X|,  1 "  Xf+4,  1 
Xi,  2"  Xl+4,  1 


_1 _ 1 

i  m  m 


fa  O' 
0  p 


xi,  r  xhi,  i  . 

xl,  2'  xf+l,  2  '  1°  ^ 


®  0  ^ '  xi.  r  xi+2.  i 


Xi,  2"  x|+2,  2 


°\  i\  i 


+  i  j|  *• 1  <+3> 1 1+  I  ^  °  I  /  X|» 1  Xi+4- 

\°  a)  XI,2“XI+  3.2  \°  “J  [  Xl„  2  ~  X  1+4.  2  / 


1  )  •  (42) 

\  /  \  ■»'  -  ij  2  / 

As  in  Section  III,  the  x  and  /  components  of  the  motion  ar^ 

completely  decoupled,  and  may  be  considered  separately.  As  a  result  of 

symmetry,  it  is  sufficient  to  consider  the  x-component  motion  for  a  complete 

quadrant  of  the  space  of  the  propagation  vector,  Jc.  The  transformation  to 

the  normal  coordinates  of  the  unperturbed  system,  Eq.  (7),  yields  the 

equation  of  motion, 


-w2^Q1  (k)exp  i  k-  +  (a/m)  ^[2-exp  (ik^aJ-expt-ik^a^Qj  (k)  exp  ik'  R^ 


+  (P/m)  }  [2-exp  (ik  a)  -  exp  ( -ik  a)]  Q ^  (k)  exp  i  k‘  R^ 


U 

k 


=  "V  (k)  [( 1 /m)-(i /m^)]  exp  i  k' R^  [a  (2-exp  i  k^a-exp  -  ik^a) 


+  (3  (2-exp  i  k^a-exp  -  ik^a)]  . 


(43) 


i 

\ 
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Lquation  (43)  is  similar  in  form  to  Eq.  (9).  By  mathematical  manipulations 
which  are  similar  to  those  which  lead  from  Eq.  (9)  to  Eq.  (11),  Eq.  (43)  is 
transformed  into 


(44) 


where 


®k,k;  *  -rr-  ^(^-•1)expi(k,-k)'R<. 


(45) 


This  result  is  identical  to  that  obtained  by  Langer  [2]  for  the  one-dimensional 
problem.  It  is  easily  shown  that  the  entire  difference  between  the  two  problems  i s 
represented  in  the  form  of  the  function  G^^u^)  which,  for  this  problem,  is 
given  by 


(1)  2  qXul,k 

G.UV)  =  - *— 


1  + 


*1 


“ltf 


(46) 


k' 


where 


K  =  (m/m1)  -  1 


(47) 


The  remainder  of  the  problem  consists  of  the  evaluation  of  the  sum 
2 


wl,k‘ 


_  4«  y 

m  [_j 


sin  (kht/2) 


T,k' 


+  ^  -  sin-(k;a/Z) 


m  [_j  2  2 

k'  U1 ,  k'  "  U 


(48) 
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The  two  sums  which  appear  cn  the  right-hand  side  of  Eq.  (48)  have  been 
evaluated  in  the  Appendix  I  and  have  been  used  in  Section  III,  for  the  case 
where  Eq.  (34)  holds.  Substituting  Eqs.  (33)  and  (35)  into  Eq.  (48),  and 
inserting  the  result  into  Eq.  (46),  one  obtains  the  result 

G[U(w2)  =  q\(02k{l+K[l  +  (i£j2m/(aP)1/2]}  .  (49) 

With  this  result,  it  is  easily  shown  that  the  frequencies  of  propagation, 
in  the  elastic  limit,  are  given  by  the  expression 

U2  =  W2  k  (1  +  qK  )  ,  (50) 

where 


K  -  1  -  (m'/m)  • 

The  inverse  lifetimes  are  given  by 

2  2 

(1)  qX  mu  ls  k 
^  2(l+\)2(ap)1/2 


(51) 


(52) 


Before  leaving  this  problem,  we  should  point  out  that  another  aspect 
of  it  has  been  treated  by  Mahanty,  Maradudin,  and  Weiss  [3].  They  considered 
the  local  modes  which  arise  when  m'  <  m.  Many  of  the  mathematical 
difficulties,  which  these  authors  surmounted  quite  elegantly,  have  been 
avoided  in  this  work  because  attention  was  restricted  to  the  elastic  spectrum. 


V.  DISCUSSION  -  THE  VIRTUAL  CRYSTAL  APPROXIMATION 

One  frequently  finds  statements  in  the  literature  that  the  so-called 
"virtual  crystal  approximation"  must  hold  for  the  elastic  spectrum  of  a  crystal 
containing  a  small  number  of  impurities,  or  a  small  degree  of  disorder.  In 
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this  approximation,  the  masses  are  all  taken  as  the  average  mass,  and  the 
spring  constants  are  all  taken  as  the  average  spring  constant  The  justification 
given  for  this  approximation  is  usually  that  the  wavelengths  of  the  modes  are 
very  long  compared  to  the  range  of  the  local  perturbations,  so  that  these 
perturbations  are  just  averaged.  Other  authors  have  carried  the  approximation 
further,  and  have  concluded  that  the  failure  of  the  virtual  crystal  approximation 
to  explain  the  vibrational  frequencies  near  the  Brillouin  zone  surface,  of 
dilute  alloys,  is  evidence  for  the  existence  of  short-range  order  in  these 
alloys.  In  this  section,  we  examine  this  approximation  with  the  aid  of  the 
solutions  obtained  in  the  preceding  sections. 

The  two  one-dimensional  solutions  obtained  previously  [1,2]  yield 
virtual  crystal  approximation-type  results  to  first  order  in  q  Langer  [l] 
finds  that  the  elastic  spectrum  of  a  system  with  is  topic  impurities  is 
completely  described  by  an  average  mass  given  by 

m  =  (l-q)m  +  qm1  .  (53) 

In  (P-B),  the  authors  find  that  if  is  the  inverse  of  the  spring  constant 
which  must  be  averaged.  That  is,  taking  y  as  the  nearest-neighbor  spring 
constant  of  the  unperturbed  system,  and  7'  as  the  spring  constant  which 
replaces  it  at  randomly  positioned  sites,  the  elastic  spectrum  is  completely 
described  by  an  effective  spring  constant,  7  ,  which  is  given  by  the  expression 

1/7  =  [<1  -q)/7]  +  (q/V)  •  (54) 

Henceforth,  we  shall  speak  of  this  equation  as  the  virtual  crystal  approximation 
for  spring  impurities. 
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For  the  isotopic  substitution  in  two  dimensions,  the  situation  is  quite 
similar.  Equal  .n  (50)  can  be  rewritten  as 

if  =  m  10  2  k  [m(l-q)  +  qm1]  1  (55) 

2 

to  first  order  in  q  .  Since  k  is  independent  of  mas 3,  this  is  the  result 

expected  from  the  virtual  crystal  approximation. 

There  is  good  reason  to  believe  that  the  approximation  must  held  for 
the  elastic  spectrum  of  a  system  of  any  dimensionality  and  any  concentration 
where  the  only  substitutions  are  isotopic  and  the  system  is  homogeneous  on  a 
macroscopic  scale.  The  reasoning  goes  as  follows.  The  masses  enter  the 
dynamical  equation  through  acceleration  terms  only.  Thus,  they  do  not 
contribute  at  all  to  the  static  elastic  constants,  which  depend  only  on  the 
interactions.  The  only  way  that  the  masses  enter  into  the  classical  expressions 
for  the  elastic  spectrum  is  through  the  density,  since  the  elastic  spectrum 
is,  of  course,  entirely  determined,  for  a  homogeneous  solid,  by  a  knowledge 
of  the  elastic  constants  and  the  density.  Thus,  the  only  condition  for  the 
applicability  of  the  virtual  crystal  approximation  to  the  elastic  spectrum  of 
an  isotope-substituted  system  is  that  it  be  homogeneous.  That  is,  one  must 
be  able  to  divide  the  crystal  into  volume  elements  whose  dimensions  are  small 
compared  to  the  wavelength  of  the  elastic  modes,  all  of  which  have  the  same 
density  associated  with  them.  In  this  case,  the  virtual  crystal  approximation 
is  just  a  means  of  calculating  that  density. 
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The  argument  just  presented  leads  one  to  suspect  that  it  is  in  the  case 
where  the  spring  constants  are  changed  that  one  expects  a  breakdown  of  the 
approximation.  This  is  the  case.  For  comparison,  we  may  rewrite  Eq.  (54) 
as 

Y  =  7  [1  +  q^T/(T  +  *'  7)]  •  (56) 

Comparison  of  Eq.  (38),  for  the  perturbed  system,  with  Eq.  (4)  for  the 
unperturbed  system,  shows  that  one  can  define  effective  force  constants,  a 
and  (3  ,  by  the  relations 

—  _  ,  .  qAa 

(X  -  Qt  1  .  ^  \  l  J -} 

a  +  Aa{Z/n)  tan  (  /jB ) 

and 


(57) 


j5-  =  p 


l  +  _  qAp _ 

—  p  +  Ap  (2/tt)  tan11  (p/a)T/^ 


(58) 


These  results  indicate  immediately  that  a  knowledge  of  q  ,  a  and  a'  is  not 
sufficient  to  define  a  .  The  relationship  between  a  and  p  also  enters  into 
the  expression  for  a.  The  same  is  true  for  p  . 

The  reasons  for  this  dependence  can  be  visualized  easily.  Consider 
Fig.  1,  which  portrays  the  displacements  of  atoms  in  the  region  of  an  impurity 
spring  from  equilibrium  while  participating  in  a  longitudinal  mode.  Here,  we 
have  taken  «'  >  a  .  The  intersections  of  the  dashed  lines  represent  the 
equilibrium  positions.  The  straight  vertical  lines  represent  the  lines  of 
constant  phase  associated  with  the  mode  of  the  unperturbed  system.  Atoms  1 
and  2  are  connected  bv  a  spring  which  is  stronger,  and  which  prevents  them 
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from  separating  by  the  distance  associated  with  the  unpe^tv  ibed  mode.  As 
a  result,  there  is  a  shear  component  in  the  motion  of,  say,  atom  J  with 
respect  to  atoms  3  and  5,  and  the  she&r  constant  (3  enters  into  the  frequency 
of  a  mode  whose  frequency  is  determined  entirely  by  a  in  the  unperturbed 
system.  Some  of  the  relationships  of  this  picture  to  the  mathematical  solution, 
Eq.  (57),  are  easily  seen.  We  note,  first  of  all,  that  this  is  a  calculation  to 
first  order  in  q  .  Thus,  we  can  neglect  other  "wrong"  springs.  As  a  result, 
the  picture  indicates  that  there  is  no  shear  motion  of  atom  1  with  respect  to 
atom  2.  Thus,  |3'  does  not  enter  into  the  expression  for  a  .  Similarly,  we 
would  expect,  from  this  picture,  that  the  virtual  crystal  approximation  would 
hold  for  a  if  {3  =  0.  This  is  the  case,  since  (Z/u)  tan  *(oo)  =  1,  making 
Eq.  (57)  identical  to  Eq.  (56).  Setting  (3  =  0  also  has  the  effect,  of  course, 
of  making  the  lattice  unstable.  Finally,  we  may  look  at  the  problem  as  the 
picture  does,  with  a  classical  perturbation  theory  viewpoint,  treating  (A a) 
as  the  perturbation  on  this  mode.  In  first  order,  the  energy  shift  is  obtained 
by  letting  the  system  undergo  its  .nperturbed  motion.  In  this  order,  j3  should 
not  appear.  It  is  only  in  second  order,  when  the  motion  has  been  changed  to 
first  order,  that  the  effect  of  (3  should  ente~  in.  It  is  readily  seen  from 
Eq.  (57)  that  |3  only  enters  into  terms  of  second  and  higher  order  in  A  a  . 

Equations  (£7)  and  (58)  do  show  that  the  two  virtual  crystal  approximations 
which  may  be  used  to  obtain  average  spring  constants  form  upper  and  lower 
bounds  for  the  actual  a  and  (3  ,  for  the  model  discussed  here.  The 
tan  ^  (a/P)^  ^  and  tan  *(p/<*)*^ 


in  these  equations  have  upper  and  lower 
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bounds  of  0  and  n  /2 


_  | 

For  tan 


0  ,  these  equations  become 


1  /at  [ s  i  -q> fn>\  f  [q/V]  (59) 

and 

T  =  (l-q)P  +  qP'  .  (60) 

For  tan  ^  =  77/2,  of  course,  the  situations  are  reversed.  Equation  (59) 

corresponds  to  the  l/y  averaging  which  was  obtained  in  one  dimension,  while 

Eq.  (60)  is  the  straightforward  averaging  of  the  spring  constants.  Any  inter- 

- 1  1/2 

mediate  values  of  tan  A  (f i/o)  '  will  yield  results  whic*>  lie  between  the  two 
simple  forms  of  averaging.  Note  that  this  discussion  has  been  restricted  to 
the  cases  where  Aa/«>  -land  A(3/{3  >  1.  These  conditions  must  be  satisfied 

if  each  atom  has  its  motion  centered  bout  a  position  of  stable  equilibrium 
in  the  static  crystal. 

Paul  [4]  has  obtained  the  identical  upper  and  lower  bounds  for  the 
elastic  constants  of  isotropic  multiphase  materials.  He  shows  that  such  bounds 
hold  in  general  in  the  situation  where  .  .  the  constituents  are  distinct  and 
capable  of  separation  by  purely  mechanical  means  (e.  g.  ,  no*  solid  solutions)". 

In  the  language  of  lattice  dynamics,  Paul  is  restricting  his  work  to  situations 
in  which  the  force  constants  between  atoms  belonging  to  the  interiors  of 
different  constituents  are  zero.  That  is,  the  elastic  constants  associated 
with  each  of  the  constituents  completely  describe  the  interactions.  This 
situation  is  quite  similar  to  that  discussed  here.  In  the  determination  of  the 
e!as“ic  cons^an's,  it  is  only  the  pairwise  interactions,  and  not  the  masses, 
which  are  significant  Thus,  since  nearest  neighbor  interactions  only  are 
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considered,  we  may  consider  each  impurity  spring  as  an  independent  constituent, 
and  the  results  are  consistent  with  Paul's  bounds.  For  some  reason  that  we 
do  not  understand,  the  fact  that  the  system  is  not  isotropic  does  not  seem  to 
be  important  in  this  case.  This  analysis  raises  the  possibility,  however, 
that  important  and  interesting  effects  may  appear  when  longer  range  inter¬ 
actions  appear. 

Figure  2  shows  the  deviations  of  a  and  ]3  ,  as  given  by  Eqs.  (57)  and 
(58),  from  the  virtual  crystal  approximation.  In  these  graphs,  we  have  taken 
ar/p  =  9.  This  ratio  was  obtained  from  the  elastic  constants  of  j3-brass,  as 
measured  by  McManus  [5]  by  setting  a/\ 3  =  /[(C^  -  C^)/2].  As 

expected,  a  lies  close  to  the  curve  defined  by  Eq.  (59),  while  (5"  is  given, 
approximately,  by  Eq.  (60).  A  choice  of  ar/6  ~  1,  corresponding  to  a 
strongly  covalent  crystal,  wculd  lead  to  lines  for  a  and  which  have 
slopes,  at  q  =  0  and  1,  that  are  approximately  equal  to  the  average  of  the 
slopes  associated  with  the  two  virtual  crystal  approximations. 

Since  a  and  ^  also  determine,  through  Eq.  (41),  the  spectral 
distribution  function  for  small  ,  the  remarks  about  a  and  ^  apply  to 
it  as  well. 

Finally,  it  is  important  to  note  that  the  random  spring  calculation 
performed  in  this  work  applied  specifically  to  the  order -disorder  problem. 
There  exists  another  interesting  related  situation  which  remains  to  be 
treated.  It  is  the  case  where  impurity  atoms,  which  are  substituted  into 
the  two-dimensional  lattice, have  approximately  the  same  mass  as  the  solvent 
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atoms,  but  change  the  spring  constants.  In  this  case,  the  four  sets  of  spring 
constants  associated  with  the  interaction  of  the  impurity  atom  with  each  of  its 
neighbors  must  appear  together  in  the  averaging.  Work  is  underway  on  this 
problem,  and  it  is  expected  that  the  results  will  soon  be  presented.  The 
results  obtained  here  indicate  that  it  is  highly  likely  that  the  virtual  crystal 
approximations  are  inadequate  for  the  description  of  the  modifications  of  the 
elastic  spectrum  of  such  a  system.  If  this  is  the  case,  any  attempt  to  assign 
observed  breakdowns  of  the  approximation  to  ordering  or  clustering,  without 
corroborative  diffraction  evidence,  would  appear  unjustified.  This  is 
particularly  true  when  the  modes  are  not  a  portion  of  the  elastic  spectrum, 
but  are  associated  with  the  Brillouin  zone  surface. 

VI.  ACKNOWLEDGMENTS 

This  work  was  supported  in  part  by  the  National  Science  Foundation 
under  Grant  GK-222,  by  the  Advanced  Research  Projects  Agency  under 
Contract  SD-88,  and  by  the  Division  of  Engineering  and  Applied  Physics, 
Harvard  University,  Cambridge,  Massachusetts. 

The  authors  are  indebted  to  U.  F.  Kocks  and  B,  Budiansky,  who 
pointed  out  the  existence  of  the  work  of  B.  Paul  [4]. 


ARPA-18 


2-28 


APPENDIX  I 

EVALUATION  OF  THE  INTEGRAL  OF  EQUATION  (32) 
The  value  of  the  integral 


jr/a  nj a 


I  = 


sin2(k  a/2)  dk  dk 
_  x  /  x  y 


11m 

-ff/a  -*i/a  (4a/m)sin2(k  a/2)+(4(3/m)ain2(k  a/2)-u2-ie 

x  y 


(A  - 1 ) 


is  required  in  the  limits  y2«4cr/m  and  u2«  4{3/m.  We  begin  by  fixing 

2 

ky  and  evaluating  the  integral  over  k^.  Consider  the  case  w  >(4|3/m) 
sin2(k  a/2).  Let 

y 

,2  2 
1  —  t  t 


u'  =  u4  -  (4^/m)  sin2  (k^a/2)  , 


(A  -2) 


and 


*nn4(k  a/2)  dk 


•  2, 

.  a  r  ,In  (  x  _ 

1  _  (4a/m)  sin2(kxa/2)  -  w'2-  ie 


Letting 


z  =  exp  ik  a  , 


u"2  =  (m/ar)u*2  , 


and 

e"  =  (m/or)  e  , 
Eq.  (A-3)  becomes 


m 


'1  4aai 


5 


( z  - 1 )  *  dz _ 

z  [z  -(2-0)”  -ie")  z  +  :  j 


(A -3) 


(A -4) 
(A -5) 


(A-6) 


(A-7) 
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where  the  contour  of  integration  is  the  unit  circle.  The  denominator  of  the 
integrand  has  three  roots, 


z0  "  0  » 


1  ..2  1.  1  1 


( A  -  8 ) 

u"2~  jU’1)2]1/2 

(A-9) 

,»2-‘  ic")2]1/2 

(A- 10) 

z2  =  {l~2u "  -y1*")  "  1  [l-O-  ’?  i€'">  ]  (A-10) 

Since  u"2  <1,  z^  is  inside  and  z^  is  outside  the  unit  circle.  Thus,  we  obtain 
the  solutions 


I  -  IB1  J  1  +  (zi-1) 

1  Zaa.  z  i  z  2  2  i  (  z  |  -  2  £ 


(A-  1 1 ) 


Evaluating  Eq.  (A- 11)  in  terms  of  the  variables  of  Eq.  (A-l),  we  obtain 
m7r  f  i(p/a)l/2[(mu>Z/4P)-Sin2(kya/2)]l/Z 
1  Jaa  [l-(mU2/4or)  +  (P/or)  sin2  (k^a/2) /2  J 


(A-12) 


z  z 

In  a  similar  manner,  we  consider  the  case  u  <  (4|3/m)  sin  (kya/2). 


Letting  u'2  =  (4(3/m)  sin2  (kya/2)~(j2  , 
we  obtain  the  integral 


(A-l  3) 


sin2(k  a/2)  dk 

L  =  1  _ 2L.— _ * _ 

2  (4a/m)  sin2(kxa/2)+u' 2 


(p /or)  *  ^2[sin2(k  a/2)-(cj2m/4|3)]^2 


[l+((3/a)sin  (kya/2)-(u  m^o')]  ' “  . 


(A- 14) 


Substitution  of  Eqs.  (A-12)  and  (A-14)  into  Eq.  (A-l)  yields  the  expression 
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I  = 


77/a 

am  i  dk 


6 77 or  J  Vr 

-77/a 

iam  <2A>«-W/4m1/2  (pA)i/2[(^m/4p,.sin2(k  /2)]i/2dki 

+  j - x- 


am 

lira 


0  [l-{mw2/4Qr)+(p/o')sin2(kya/2)] 


TJT 


it  /a 


(p/a)1//2[sin2(k  a/ZJ-tw^/^)]1/4'  dk 


_ y  f  y 

-  _  2  ,  /?  [l+(P/af)sin2(k  a/2)-(u)2m/4a)]1^2 

(2/a)sin~  (mu  /4|3)  '  ’ 


(A- 1 5) 

It  should  be  noted  at  this  time  that  tha  integrands  are  all  real.  Thus, 
the  first  and  third  terms  of  the  right-hand  side  of  Eq.  (A- 15)  contribute  the 
real  part  of  I,  while  the  second  term  yields  the  imaginary  part. 

For  convenience,  we  denote: 


TT/a 


I-  =  (am/STra)  \  dk  =  m/4a  , 

-f/.  y 


(A- 1 6) 


iam 


4i3"( ~  } 


(2/a)8in-1(mw2/4p)i//2  ?  ?  1 /? 

j^_^l/2  f  [(mu  /4j3)-sin~(kya/2)]  /  dky 

0  [  1  -(mu'/4a)+(^ /a)sm^{k  a/2)j  ^2 


(A-17) 


i  -  (  JL*1/2 

l5  "  4ira  {  a  ’ 


77/a 

I 


[sin2(kya/2)-(mu2/4p)]1^2  dky 

2  l  A  \  .  \  !  2/,  /->  V  l  1  j 2 


.i  2  i  /?  [ln(mw  /4or)+(p/ar)sin  (k  a/2)] 
(2/a)sin~  (mio  /4p)  '  y 


(A- 18) 
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t  =  sin  (ky.a/2)  , 


6  =  (mu  /4|3)  «  1  , 


6'  =  (mu  /4 a)  «  1  . 


(A  - 19) 


Then  Eq.  (A- 17)  can  be  written  as 

,  (6)‘/2 


T  -  ima 
4  8 77a 


{JL)1/2  7  (S-t2)1/2^ 

*  *  ri-6'  +  (8/or)t 


(8/a)t"]1/"[l.-t"]J 


(A-20) 


It^«  easily  seen  that  t  «  l  over  the  range  of  integration.  Thus,  we 
may  approximate  the  entire  denominator  of  the  integrand  of  Eq.  (A-20)  by  unity 
to  obtain  an  expression  which  is  good  to  only  first  order  in  5  and  5  '.  In 
this  case  we  may  write 


0 

im2u2  ,  (3  xl/2 
-  \*>~A —  V-*r) 


(A -21) 


Substitution  of  Eq.  (A- 19)  into  Eq  (A- 18)  yields  the  expression 


5  “  2 f 


""  *  [l-t 


u2^)1/2  dt 


[  1  ‘-6'  +  (|3/a)  t  ] 


(A-22) 


In  this  case,  the  approximations  are  less  straightforward.  It  is 
dear  that  5  ’  may  be  set  equal  to  zero,  since  it  is  always  very  much  less 
than  the  other  terms  contained  in  the  radical  in  which  it  appears.  Further 


ARPA-18 


2-32 


simplification  of  the  integrand  is  hindered  by  the  fact  that  6  /t^  ^  1  near  the 

2  1  /2 

lower  limit  of  integration.  Nevertheless,  we  may  write  [l-(6/t  )]  '  ^  1  in 

2  i  !z 

the  integrand  for  the  reasons  which  follow,  The  factor  [1  I  ((3/a) t  ]  of  the 
denominator  is  a  relatively  slowly  varying  function.  As  a  result,  the  remainder 
of  the  integrand  is  sharply  peaked  near  the  upper  limit  of  integration.  In  this 
region,  it  lite  reasonable  to  replace  the  radical  of  the  numerator  by  unity. 

Hence,  if  contained  the  entire  contribution  to  the  real  part  of  I  ,  it  would 

be  sufficient  to  make  this  replacement.  Unfortunately,  however,  a  major 
portion  of  Ij.  is  cancelled  by  I^.  Thus,  it  is  necessary  to  justify  thi3  replace¬ 
ment  a  bit  more  carefully.  We  do  this  by  calculating  an  upper  and  lower  bound 
to  the  integral.  It  is  then  shown  that  the  difference  between  the  two  is  negligible 
compared  to  the  real  part  of  I.  To  compute  these  bounds,  we  note  that  the 
integrand  is  positive  and  real  over  the  range  of  integration.  In  addition,  over 
this  range,  [  1  -{5/i^)j  <  [1  -(6 ,/t^)]*^  <  1 .  Thus,  lower  and  upper  bounds  to  the 
integral  can  be  obtained  by  using  the  upper  and  lower  bounds  of  the  radical 
in  the  integrand.  To  obtain  the  values  of  the  associated  integrals,  we  need 


(A-23) 
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APPENDIX  II 


EVALUATION  OF  THE  SPECTRAL  DISTRIBUTION  FUNCTION 
FOR  THE  IMPURITY  SPRING  PROBLEM 
As  indicated  in  Eq.  (17),  the  spectral  distribution  function  for  this 
problem  is  given  by 

^(u)=T  N-^oo  TT  taTr  D(1)(uZ+ie).  (A-26) 

e-0 

Using  Eq.  (26)  for  D^  and  Eq.  (36)  for  G^,  Eq.  (A-26)  becomes, 

2 

for  small  u  . 


g  ( w  )  =  (4u/Nir)T  [(4Q-/m)sin2(kxa/2)+(4p/m)sin2(kya/2)-w  2-ie]~ 1 
€~*"0  k 


=  (a2u/ 7T3)  Im 


j  j  [(4ar/m)sin2(kxa/2)+(4p'/m)sin':(ky.a//2)-u':-ie]“Adkcdky. 
-IT /a.  -7T/a  (A -27) 


•  2, 


2  .  -,-1 


Here,  a  and  (T  are  defined  by  Eqs.  (57)  and  (58),  respectively. 

Let  J(u2)  be  the  integral  ox  Eq.  (A =27).  It  can  be  written  as 

I  =  2(1,  +  I2)  , 


(A-28) 


where 


(2/a)sin”*(mu2/4p)*/2 


I 


dk 


IT/a 


[(4a/m)sin4'(k  a/2)+(4p/m)sinfc(k  a/2)-u4-ie]" *  dk  ,  (A-29) 

x  y  x 


.2  .  n-l 


■ir/a 


x 
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[2  - 


[(4a'/m)sin^{kxa/2)+(4^’/m)sin^(kya/2)-u)^-i€] dkx  .  (A-30) 


The  real  part  of  the  denominator  of  the  integrand  of  is  always  positive.  As 
a  result,  the  imaginary  part  of  the  integral  is  of  first  and  higher  order  in  e  , 
and  vanishes  when  the  e  limit  of  Eq.  (A- 27)  is  taken. 


f2  =  mu2/4jT  , 

U)"2  -  (m/a  )  [w2-(4^/m)  sin2  (k^a/2)]  , 


-  f  1 


a  -i/ a  4  pin  <kxa/2)*w"  -i€ 


(A -31) 


(A-32) 


(A-33) 


Then 


(2/a)sin_1  £ 


l,  dk  . 

3  y 


(A-34) 


Substitution  of  Eq.  (A-4)  into  Eq.  (A-33)  yields 


im  f  dz 

—  i  7  2  > 

aa  ^  z  -  (2-u"  -ie)z  +  1 


(A -3  5) 


where  the  contour  is  the  unit  circle. 


r 
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The  denominator  of  the  integrand  of  Eq.  (A-35)  is  similar  to  that 
of  Eq.  (A-7).  It  has  roots,  z^  and  z^  ,  given  by  Eqs.  (A-9)  and  (A-10), 
respectively.  As  before,  z^  lies  inside  and  lies  outside  the  unit  circle. 
Thus, 


I3  =  (mi/aa)  [2iri/(z1  -  z2)] 

=  (Zvi/a.)  {[u2-(4j3/m)  sin2  (k  a/2)] 

5  .  .  1  /■> 

(A-36) 

After  some  simple  manipulations  which  use  Eq.  (A-31),  substitution 
of  Eq.  (A-36)  into  Eq.  (A-34)  yields 
(2/a)sin'1f 


x  [(4or/m)  +  (4p/m)  sin2  (k^a/2)-u2] }  *^2 


I,  = 


unir 


1  "  2a(«F)l/2 


I 

0 


dk 


[?2-sin2(kya/2)]^2[l-(mw  /4o-)+(p/a)sin2(kya/2)]1  ^2 


(A -37) 


To  lowest  order  in  u  ,  v/e  have 

sir 

I 


(2/a)sin_1  | 


im7T 


dk 


2a(o-p):/2 


[?2-  sin2(kya/2)]^2 


im7T 


2a  2(cr^)1//2 


(A -38) 


Substitution  of  Eqs.  (A-28)  and  (A-38)  into  Eq.  (A  -27)  yields  the 
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FIGURE  CAPTIONS 


Fig.  1  Pictorial  representation  of  the  displacements  of  atoms  in  the 
vicinity  of  an  impurity  spring,  while  taking  part  in  a  (10) 
longitudinal  mode.  The  dashed  lines  represent  the  equilibrium 
positions.  The  impurity  spring,  with  Ot, 1  >  a  ,  is  between  the 
atoms  denoted  1  and  2.  Shear  along  the  lines  3-1-5  and  4-2-6 
introduces  the  factor  3  into  the  expression  for  the  effective  C.  . 

Fig.  2  Comparison  of  the  effective  force  constants,  a  and  fT,  as  a 

function  of  concentration,  q,  with  virtual  crystal  approximation 
predictions.  Here,  A  Cl/  &  and  A  3/ 3  have  been  taken  as  -0.  3, 
while  a/  3  =  9.  The  dashed  line  represents  the  linear  averaging 
of  the  force  constant,  while  the  solid  curve  shows  the  averaging 
of  the  inverse  of  the  force  constant. 


FIGURE  1 


